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Distributed source identification for wave equations : an 

observer-based approach* 

Marianne Chapouly and Mazyar Mirrahimi ^ 

Abstract 



In this paper, we consider the ID wave equation where the spatial domain is a bounded 

interval. Assuming the initial conditions to be known, we are here interested in identifying 

^H an unknown source term, while we take the Neumann derivative of the solution on one of 

" the boundaries as the measurement output. Applying a back-and- forth iterative scheme 

' and constructing well-chosen observers, we retrieve the source term from the measurement 

0^ output in the minimal observation time. 
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'5 1 Introduction 

a 

H Blum and Auroux [l] proposed a new inversion algorithm for identifying the initial state of an 

'— ' observable system, based on the application of back-and-forth observers. Noting that we have 

_. •• only access to a measurement output on a fixed time interval (0, T), the idea consists in proposing 

^j^ a first asymptotic observer for the system that will be applied in this time interval and a second 

00 one that will be applied to the system where the direction of time is reversed. These observers 

C^ are then used iteratively to get a better estimate of the initial state after each back-and-forth 

^^ iteration. If the two observer gains are well-chosen, so that the whole back-and-forth procedure 

i induces a contracting error dynamic, one can ensure the convergence of the estimator towards 

•O the initial state. 

More recently, Ramdani-et-al [H [S] have considered, for the case of wave equations, the 
theoretical study of this problem applying techniques borrowed from semigroups theory. Here, 
we consider a similar problem to [S] for a wave equation, where the initial state is known (is zero) 
and we are rather interested in identifying an unknown source term q{x). Let also T > 0, a; G M 
and let q g -ff^(0, 1) n -ffo(0, 1). We consider the following system 

H 

"^ ( utt- u^x = q{x) cos(a;i), (t, x) G (0, T) x (0, 1), 

I w(i,0)==u(i,l) = 0, tG(0,T), 

^' I w(0,x) =ut(0,a;) = 0, xG (0, 1), 

[ 2/(t) = ".(t,0),<G(0,r), 

where (u, ut) represents the state of the system and y is the output. The term cos{ujt)q{x), where 
w is a fixed (known) frequency and q{x) is unknown, is considered to be an external force which 
varies harmonically. 

For any q G i?^(0, 1) n H^{0,1) there exists a unique solution u G C^{[0,T];H^{0,1)) n 
C2([0,T];L2(0, 1)) to ^ and u^{t,0) G H^{0,T) (see §, Remark 2]). It is moreover well-known 
from the work [6] of Yamamoto that this problem is well-posed in the sense that one can retrieve 
the source term q{x) from the measurement of y on the time interval (0, T) if T is large enough. 
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Our aim here is to propose well-chosen observers which allow, using a back-and-forth procedure, 
to retrieve q from y in the minimal observation time. More precisely, we prove the following 
result 

Claim 1 We can construct efficient observers for which the back-and-forth algorithm is con- 
vergent and which allow, using the measurement output y(t) over the time-interval (0,2), to 
reconstruct the unknown source term q{x). 

Remark 2 Let us point out that since the spatial domain is given by (0, 1), the minimum ob- 
servability time is given by T ^ 2. 

Note that, whenever the whole initial state (w(0, .),ut(0, .),(?) is unknown, system (fTj) is not 
observable. In order to realize this, one can consider the simpler case where the source term q is 
given by only the two first modes of the wave equation, q{x) = qi sin(7ra;) + 172 sin(27ra;), where 
qi and (72 are the unknown scalars. In this case, M becomes equivalent to two independent os- 
cillators with different frequencies and with two unknown source terms qi cos(a;i) and 92 cos(a;t). 
Moreover, the output is given by a linear combination of the position of the oscillators. This six 
dimensional system with one output is not observable. However, if we know the initial state of 
the oscillators, the two parameters qi and q2 become identifiable. 

The back-and-forth estimator allows us to take into account this knowledge of the initial state 
(u(0, .), Ut(0, .)) = (0, 0) of the wave equation. On the contrary, if we had only used a forward 
observer, we would have lost the information on the initial state of the system and therefore there 
would have been no reason for the observer to converge to the real parameters. 
We prove the convergence of the algorithm using Lyapunov techniques and LaSalle's principle. 
One of the main difficulties comes from the fact that the precompactness of the trajectories is 
not ensured since we deal with an infinite dimensional system and we also have to use sharp 
mathematical estimates. 

The paper is organized as follows. 
In Section 2, we prove the equivalence between M and another system which consists in a system 
composed of a wave equation without any source term and an oscillator in cascade. Here the 
unknown term to retrieve is the initial condition. Section 3 is devoted to the proof of Claim [l] 
and is divided in two parts. In the first one we study the well-posedness of the proposed observer 
and in the second one we prove its convergence through the back-and-forth procedure. In Section 
4, we provide some numerical simulations illustrating the efficiency of the method. Finally in 
Section 5, we propose and analyze the extension of the method to the case of wave equations 
with iV-dimensional spatial domain. 

2 An equivalent estimation problem 

We begin with introducing the following system, 

wu-w.,^,^0, {t,x)£ (0,T)x (0,1), 

(2) { u;(i,o) = w(t,i) = o, te(o,r), 

w{Q, x) = q{x), wt{0, x) = 0, X e (0, 1), 

zi{t) = Z2it),t€iO,T), 
. Z2it) = -uj^zi{t) + w,it,0),te{0,T), 
^' ^ zi(0)=y(0),z2(0) = y(0), 

r(t) = zi(i),te(o,T), 

where {w,Wt, zi, Z2) represents the state of the system and Y is the output. Here, the fact that 
y{0) and y{0) are well-defined is a consequence of the definition of y (see M) and the regularity 
property ( |10[ ). System (p])-® is nothing but a homogeneous wave equation and an oscillator in 
cascade. The unknown q is now the initial datum of the system. The aim of this section is to 
prove both following results 



Proposition 3 There exists a unique solution 

{w,wt,z^,Z2) e C([0,r];iJi(0,l)) X C([0,r];L2(o,l)) x H^{Q,T) x H\Q,T) 
to ([2| . Moreover, w satisfies the following hidden property 

(4) w^{t,0)eL^{0,T). 
and 

Proposition 4 For any T > 0, for any q E H^{0, 1), 

y^Y inij2(0,T), 
where y denotes the output of system M and Y denotes the output of system (pi) . 

We start by proving the Proposition [s] which ensures the weU-posedness of the system (l2])-(|3|. 
Proof of Proposition [sl It follows from [3 Lemma 3.6 page 39] that there exists a unique 
we C([0,T];Hi(0,l))nCH[0,T];L2(0,l)) solution of 

f wtt~w,^^0, {t,x)e (0,T)x (0,1), 

(5) <^ u;(i,0) = u;(t,l)=0, te(0,r), 

[ w{0,x) ^ q{x), WtiO,x) = 0, X e (0,1). 

From the hidden regularity property (see [^ Theoreme 4.1 page 44]), 

(6) w^{t,0)eL^{0,T). 

Defining the matrix yl := 2 n ) ' ^^'^ unique solution of the equation ([3]) is given by. 

Using Q, one easily gets 

\{zi,z2)\mo,Tr < c(i + \w.j;{t, o)|l2(o,t)), 

where C denotes a positive constant which only depends on uj and T. Using first and second 
equations in ([3]), together with ^ and the last inequality we get 

(8) zi and Z2 e H\0,T), 

which implies, from the first equation in (Is]), zi G H^{0,T). This ends the proof of Proposition 

[3 ■ 

Now we start with the proof of the proposition [4| which ensures the equivalence (in the sense 
of the measurement output) between the equation TTl) and the system (l2])-(|3|. 
Proof of Proposition |4} We first prove the following lemma 

Lemma 5 The solution u to M satisfies 

(9) u e C([0,T]; ij3(0, 1) n i/i(0, 1)), "t € C{[0,T];H^{0, 1) n H^{0, 1)), 
and 

(10) u,(i,0)ei/2(0,T). 



Proof of Lemma [sl We proceed as in [6l Appendix page 495] . We recall that from [6] , we 
already know that u e C\[0,T]; H^ {0,1)) n C^{[0,T]; L^ (0,1)) and u^it^O) e H^{Q,T) (see in 
particular [6i (1.3) page 183]). We consider now the following system, 



(11) 



Vtt - Vxx = -UJ^qix) cos{ujt), (i, x) G (0, T) x (0, 1), 
v{t,0) =w(i,l) = 0,t€ (0,T), 
vi0,x) = qix),vti0,x) = 0,x£ (0,1). 



Since q G H^{0, 1) n H^{0, 1) it follows from [3, Lemma 3.6 page 39] that 

(12) z;eC([0,T];iji(0,l))nCi([0,T];L2(0,l)), 
and 

(13) w,(t,0)eL2(0,T). 
We can set 



u(t, a;) := 



^0 



v{T,x)dTds, te (o,r), xe (0,1), 



so that u e C2([0, T]; i?o (0, 1)) n C^{[0, T]-L^{0, 1)) and u satisfies ^. By uniqueness of weak 
solution to (IT]), we get 



Consequently 
(14) 



and from (13), 



u(t, x) — u{t, x) and Utf{t, x) = v{t, x), t G (0, T), x S (0, 1). 

ueC\[0,T];H^{0,l))nC'{[0,T];L'{0,lj), 

u,{t,Q)^H\Q,T). 



Now, using the fact that u satisfies the equation (fTl) and that Uu G C([0,T]; i/o(0, 1)) (see (14)) 
and that q S -^^(0, 1) n iJ,J(0, 1), we obtain 



and thus 



w,, eC([0,r];iJoi(0,l)), 



ueC([0,T];ij3nFi(0,l)). 



Using one more time the fact that u satisfies (fTl), that q e H^{0, 1) H i?o(0, 1) and that Uu G 
CH[0,T];L2(o, 1)) (see ^), it follows 



(15) 



".. eCi([o,r];i2(o,i)). 



Consequently, it follows from (14) and (^15| that 



«eC7i([o,r];ij2(o,i)ni/i(o,i)). 

This ends the proof of Lemma [5J ■ 

Since u G C{[0,T]; H^{0,1) n H^iO,l)) n C'^{[0,T]; H^{0,1) n H^{0,1)) and g £ H^{0,l)n 
iJc;(0,l), we can define W € C{[0,T], H^{0,1)) nC\[0,T]; L^{0,1)) as 

(16) W :— u u + Utt = ijJ u + Uxx + cos{ujt)q{x). 

It is easy to see that W satisfies 



(17) 



Wtt -Wxx ^ 0, {t,x) e {0,T) X {0,1), 

wit,o) = w{t,i)^o,te (o,T), 

WiO, x) = q{x), WtiO, x)^0,xe (0, 1), 



and thus, by uniqueness of the weak solution of the previous system, 

w = W mCi[0,T],H^iO,l))nC'i[0,T];L\0,l)) 

and 

w,it,0)^W,{t,0)inL^{0,T). 

In other words, 

(18) w^uj^u + u,, + cos{u;t)qix) in C([0, T],H^{0, 1)) n ^^([0, T];L^{0, 1)), 



and one can compute from (16) (computations making sense from (lOl), 

w^{t, 0) = uj^u^it, 0) + uuAt, 0) in ^^(o, T), 
and thus from (fTl) 

(19) w^it, 0) = uj^yit) + y{t) in ^^(o, T). 
Now, from ([s]), 

(20) Y + uj^Y = w^{t, 0) in L^{0,T). 
Finally, from (|3|, ([l9]), ([20]), 

Y + LU^Y = y + uj^y in L^{0, T), 



^^^^ i y(0)-y(0),y(0)-y(0), 

which concludes the proof of Proposition |4] by uniqueness of the solution to the ordinary differ- 
ential equation £_ + u^^ = 0, f (0) = ^(0) ==0. ■ 

Remark 6 One can note that using usual energy estimates and the fact that q G iJ^(0, 1) n 
-ffo(0, 1), one can prove more regularity for the solution {w,Wt) to ([2]). Precisely, one actually 
obtains 

w e C([0, T],H^iO, 1) n ^1(0, 1)) n C\[0,T];H^iO, 1)). 

As a consequence of Proposition^ we will now focus on system (|2])-((3]) and provide well-chosen 
observers for this sytem in order to prove Claim [11 

3 Observer design 

As we have already mentioned in the introduction, we are going to prove Claim [l] using back- 
and- forth observers. We first rewrite system (pi)-® in the following way 

w}=wMt,x)€iO,T)xiO,l), 

^ > ^ w^t,0)=w^{t,l) = 0,t£{0,T), 

w\0,x) ^q{x), u)2(0,x) =0, xe (0,1), 

Zi{t) = Z2{t),tGiO,T), 
, . , Z2it) = ~io^zi{t)+wl{t,0),teiO,T), 

^ ^ ^ zi(0)=y(0),z2(0) = y(0), 

Y{t) = z,{t),te{0,T), 

where, from above results, 

(«;\«;2,z\z2) g C([0,T];ij2(0,l)niJi(0,l)) X C([0,T];iJoi(0,l)) X i/2(0,l) X iJi(0,l). 



Let us recall that our aim is to propose an algorithm which allows to retrieve the unknown q 
from the measurement output Y on the time interval (0,T). Our idea consists in designing a 
well-chosen asymptotic observer, which ensures the decrease of a same Lyapunov function in 
the back-and-forth procedure. To this aim, we first make the previous system periodical. More 
precisely, we define (VV^^W'^, Z'^, Z"^) on ((0, +00) x (0,1))^ x (0,T)2 as the solution to the 
following periodical system 



(24) 



M/^i = W'^, (t,x) e (2fcr, (2fc+l)T) X (0,1), 
Wf = W^,, (t, x) e (2fcr, (2fc + 1)T) X (0, 1), 
Wt^ = ~W^, {t, x) e {{2k + l)r, (2fc + 2)T) X (0, 1), 
Wf = -W^^, (i, x) e {{2k + l)r, {2k + 2)T) X (0, 1), 
W^i(f,0) = T4^i(t,l) = 0, te (0,+oo), 
t W\Q,x) = q{x), W^{0,x) = 0, x e (0, 1), 



(25) 



f Zi{t) = Z2{t), t e (2fcr, {2k + i)T), 
Z2{t) == -Lj^Zi{t) + w^{t, 0), t e (2fcr, (2fc + i)r), 
^i(<) = -Z2{t), t e ((2fc + i)r, {2k + 2)T), 

Z2{t) = u^Zi{t) - W^{t,0), t e {{2k + l)T, {2k + 2)T), 
[ Zi(0) = y(0),Z2(0)=2)(0), 

for fc e N. Using previous results, one easily sees that 

{W\W^) e L^M+;H\Q,1) n Hl,{0,l)) X Lt:,{R+;Hl,{0,l)), 

and 

{z\z^)eHlM+)^HlM+)- 



Indeed, the above periodic system represents the system ( 22 )-( 23 1 on the time intervals (2fcT, {2k 



l)r), and the same system in the time-reversed manner on the intervals {{2k + 1)T, {2k + 2)T). 
In particular, one has 

J Y{t) = Y{t - 2kT), t e {2kT, {2k + 1)T), 

I Y{t) = Y{{2k + 2)T -t),te {{2k + l)T, {2k + 2)T). 

The asymptotic observer we propose is the following 



(26) 



W} ^W^, {t,x)e {2kT,{2k + l)T) X (0,1), 

W? = M^i„ {t, x) e (2fcr, {2k + 1)T) X (0, 1), 

Wt^ = -W"^^ {t, x) e ((2A: + l)r, {2k + 2)T) x (0, 1), 

w? = -wl,. {t.x) e ((2fc + i)r, {2k + 2)T) X (0, 1), 

W\t,{)) - 71(^1 W - Y{t)) + 7172(^3(0 - /o Y{s) ds),te (0,00), 

W\t,l) ^0, te (0,+oo), 

W\0,x)^W^{0,x)^0,xe{0,l), 



(27) 



Zi(t) = Z2{t) - i2{Zi{t) - Y{t)), t e {2kT, {2k + 1)T), 
Z2{t) = -uj^Zi{t) + W^{t,0), te {2kT,{2k + l)T), 
Zi{t) = -Z2{t) - i2{Zi{t) - Y{t)), t e {{2k + l)r, (2fc + 2)T), 
Z2(t)=cj2Zi(f)-Ty^i(i,0), ie ((2A: + l)r,(2fc-f 2)T), 

i3(i) = ^i(i), te (0,00), 
Zi(0) = Z2(0) = Z3(0) = 0. 



Here 71 and 72, the observer gains, are strictly positive constants to be fixed as the design param- 



eters. Before studying the well-posedness of system (26)-(27), let us introduce error equations 



The error term being defined as the difference between the observer and the observed system, 
error equations are the following 



(28) 



W} 



W^, (i, x) e (2fcr, (2fc + 1)T) X (0, 1), 



W? = Wl,, {t,x) e {2kT, {2k + 1)T) X (0, 1), 

Wt^ = -W\ {t, x) £ {{2k + 1)T, {2k + 2)T) X (0, 1), 

W? = -WL, (t, x) e {{2k + l)r, {2k + 2)T) X (0, 1), 

W\t,0) = 7iZi(i) + 7i72^3(i), t e (0,(»), 

W^{t,l) = Q, te (0,+oo), 

^-1(0, x) = -q{x), W^{0, x)=0,xe (0, 1), 



(29) 



Zi{t) = Z2{t) - 72^1 (t), t e {2kT, {2k + i)r), 
Z2{t) = -Lj^Zi{t) + W^{t, 0), t e {2kT, {2k + 1)T), 
Zi{t) = -Z2{t) - 72^1 (t), t e {{2k + 1)T, {2k + 2)r), 
Z2{t) = u^Z^{t) - Wl{t,Q), t € {{2k + 1)T, {2k + 2)T), 

i3(t) = ^i(t), ie(0,(^), 
[ Zi(0) = -2/(0), Z2(0) = -y(0), ^3(0) = 0. 



3.1 Well-posedness 

From now on, and till the end, we denote H^{0, 1) := {v G H^{0, 1) s.t. v{l) = 0}. Let us assume 
for the moment that the following proposition holds (for the proof we refer the reader to the 
Appendix) . 

Proposition 7 For any T > 0, for any 71 and 72 > 0, for any {q°,q^) G H^{0, 1) n -ffo(0, 1) x 
Hq{0,1) and any (^1,^21^3) G ^'^ j there exists a unique solution 

{v\v^) e C{[0,+oo); H\0,l) n H^{0,1)) x C([0, +00); //^(0, 1)), 

(6,6,6) e HH[0,+oo)) X l2([0,+(»)) X i/2([0,+oo)) 
to the following periodical Cauchy problem 

v} = v^, (t, x) e (2fcr, {2k + 1)T) X (0, 1), 



(30) 



«2 = vl^, {t, x) e (2fcT, (2fc + 1)T) X (0, 1), 
v] = -w2, (i, x) e ((2fc + l)r, (2fc + 2)T) X (0, 1), 
«2 = -«i„ (t, x) e ((2A: + l)r, (2fc + 2)T) x (0, 1), 
v\t, 1) = 0, «i(i, 0) = 7i6(i) + 7i726(<), t e ((2fc + l)r, {2k + 2)T), 
t «i(0,x) = gO(a;), f2(0,x) = gi(^), ^ e (0, 1), 



(31) 



ii{t) = 6(i) - 726(0, t e (2fcr, (2fc + i)t), 

6(i) = -^'6(t) + vl{t,Q), t e {2kT, {2k + 1)T), 

6(0 = -6(t) - 726(0, t e ((2fc + i)r, (2A; + 2)r), 
6(0 = ^Hi{t) - vl{t, 0), t e {{2k + i)r, (2fc + 2)r), 
6(0 = 6(0,iG(0,oo), 
6(0) = e?, 6(0) = 6", 6(0) =6°, 



where fc G N. Moreover (w"'^, w^, ^i, ^2, Cs) satisfy the following energy identities: for any t g 
(0,+oo), 



,2|f0|2 



1921^.(0,1) + 19^1^(0,1) +7iie2°P+7iC^1^?l 



(32) 
and 

(33) . .* 

+ 271720.^ j^ Ms)\' ds + J^u;'Mt)\' < CikX^^o,!) + I'z'IhMo,!) + 1??!' + 1^2°!'), 

where C^) denotes a positive constant which only depends on uj. 



One easily sees that the well-posedness of system (28)- (29 1 in 
C([0,+oo);i?2(o,l)niji(0,l))xC([0,+oo);iJ^i(0,l))xiJi([0,+oo))xL2([o,+oo))xij2([o,+cx))) 
directly follows from Proposition [71 The well-posedness of system (26)-(27l in 

Lt^M+; hHo, 1) n H^iO, l)) x L^^{R+;HX0, l)) x hUr+) x Ll,iR+) X hUr+) 



then follows (using also the well-posedness of system (24)-(25) in the same space). 

3.2 Asymptotic analysis 

The main goal of this section is to prove the following result, implying Claim [T] 

Theorem 8 For anyT>2, 

(34) lim {W\2nT,.),W^{2nT,.)) ^ {0,0) m H^iO,!) X L^{0,1), 



(35) 



lim {Zi{2nT), Z2{2nT), Z3{2nT)) = (0,0,0). 

71— >+00 



Proof of Theorem ^Let us assume that ( 34 ) and (351 do not hold. Then, there exist a positive 
constant a and a subsequence {(l>{n))n>o G N" with lim (j){n) = -l-oo such that for any n > 0, 



n— f +00 



t=24>{n)T 



(36) max {\\W\t, .)|ki(o,i), \\W^{t. .)IU^(o.i), 1^^(01, l^'WI, \Z'{t)\ 

From Proposition [7j 

[W^,W'^) eL°°(M+;i72(o,i)niJ^i(0,l)) xL°°(M+;i7i(0,l)), 

(^1,^2,^3) e H\R+) X L'{R+) X H^{R+), 
and for any time t E M.*^, 

|W^i(i,-)li.(o,i) + |W^'(^.-)li.(o,i)+7i|^2(i)P + 7i^'l^i(i)P + 27i72C.^l^ili.(o,*) 

l9lHMo.i)+7il2/(0)P+7i^'|y(0)p. 



> a. 



(37) 



Since {(W\24>{n)T, .), W'^{2(t){n)T, .)), n > 0} is bounded in H^{0, 1) x H^^O, 1), it follows from 
Rellich's theorem that there exists a subsequence of (0(n))„>o, that, for convenience, we still 
denote (0(n))„>o, and there exists (Vr°°.i, W°°.2) g ^i^q^ 1) x L2(0, 1) such that 



(38) lim {W^{2(P{n)T, .), W^{2<j){n)T, .)) = {W°°^\ W°°'^) in H^{0, 1) x ^^(o, l). 

n-4-4-00 



From (|37|, W^{2(t){n)T,0), Zi(2(/i(n)T) and Z2{2(t){n)T) are bounded, and from third and sixth 
equations in (281 so is Z^{2(j){n)T) . Thus, (up to a subsequence of (0(?t.))„>o), there exists 



(40) 



{Zf, Z^.Zf) eM? such that 

(39) hin (Zi(20(n)r),Z2(20(n)T),Z3(2(/)(n)r)) = (Z-,Zr,Z3-). 

71— J' + OO 

Let now 

{v\v\xi,X2,x^) e C([0,T];i/i(0,l)) x C([0,T]; L2(0, 1)) x L\Q,T) x L\Q,T) x ^^(o^t) 
be solution to 

i;i=«2^(i,x)e(0,r)x(0,l), 

i;2=i;i^, (t,a;)e(0,T)x(0,l), 

v^{t, 1) = 0, t;i(i, 0) = 7ixi(i) + 7i72a;3(i), t G (0,T), 

a;'i(t) = a:2(i)-72a;i(i),i€(0,r), 

i2(i) = -tj2a;i(t) + wi(t, 0), t e (0, T), 

i3(i) = 2:i(i), te(0,T), 

^1(0, a;) = VF°°'i(a:), ^^(o, x) = W^^^'^^^,)^ 3, ^ (q, 1), 

[ XI (0) = Zr, X2(0) = ^2°°, X3(0) = Z3-. 

Finally, let us define, for any n > 0, 

(z;i'",z;2.»,a:^x^,a;^) € C([0,T]; i/,i(0, 1)) x C([0, T]; ^^(o, l)) x ^^(o^r) x L\0,T) x i2(o,T) 

by 

(41) («i-"(t, .), z'''"(t, .)) := {W\2^{n)T + i, .),W\2(^{n)T + t, .)) 
and 

(42) {xm^^m^^m) ■■= (^i(2</.(n)T + <),Z2(2(/.(n)r + i),Z3(2^(n)r + i)), 

for any t £ (0,T). With such a definition, for any n > 0, {v-^'^,v^'^,Xi,X2,x'^) is solution to 

' vl'"=v^^\(t,x)e{0,T)x{0,l), 
«^"=«i'", (i,x)e(0,T)x(0,l), 

i;i'"(i, 1) = 0, v^'^it, 0) = 7ia^?(<) + 7i72a;^(t), i e (0, T), 
i?(t)=a;2"(t)-722:ni), te(0,T), 
i2"(t) = -co^xUt) + vl'^it, 0), i e (0, T), 

i^(i) = a;?(t),te(o,r), 

t;i'"(0, a;) = W^{2(j){n)T, x), x e (0, 1), 
^2,n(o^ x) = W^2(20(n)T, a;), a; e (0, 1), 

[ xUO) = ^i(20(7i)T), x^{0) = Z2(20(7i)T), x^(0) = Z3((/.(n)T). 



(43) 



From (32 1, (38l-(|40[) and (43) (more precisely by continuity of flow with respect to the initial 
state) , 

(44) lim (^;l■^^,2.") = (^;^^;2)inL-((0,^);iJl(0,l))xi-((0,^);L2(0,l)), 



(45) 



lim (a;5',4,4) = (xi,X2,X3) inL°°(0,T)= 



We now introduce the following Lyapunov function 



V(W\W^Z,,Z2,Z3m:=l 



i\W^it,x)\' + \W'{t,x)\')dx 



+ 7ic.2|Zi(i)|2+7i|Z2(t)n,te(0,r). 



Indeed, V(i) > 0, < > and using (28) and (29), one can compute, for any time t E (0,T), 
d 



(47) 



dt 



V{W\W^ Zi, Z2, Z3){t) = -7i72L^2|^i(i)|' < 0. 



The function V is positive, decreasing. Consequently there exists I > such that 
(48) hm V{W\W^,Zi,Z2,Zi){t)^l. 



On the other side, from (41), (42), (44 1 and (46), one has, for any t e (0,r) 



hrn V{W\W\W\Z,,Z2,Z3K2cl,{n)T + t) = \( ( {\vl{t,x)\'' 

(49) n^+oo 2 V Jo 

+ |w2(t, a;)|2) dx + 7i0j2|xi(t)|2 + 7i|a;2(t)|2) , 



Thus it follows from ([48]) and (|49| that for any time t e (0,r), 

V(v'^ .v"^ ,Xi,X2,X'i){t) = /. 

In other words, 1 1— ^ V{v^ ,v'^ ,xi,X2,x^){t) is constant on (0,T) and thus 



(50) 



di 



V{v\v\xi,X2,x^){t) = -7i72W^|xi(t)r = 0, i e (0,r). 



Therefore, we obtain 

(51) xi = on (0,r). 
Consequently, from the fourth and sixth lines in ( [40| , we also get 

(52) X2 = Q, X3 = Z^ on(0,r), 



where Z|° is a real constant. Then, using the third and fifth lines in (40) and, (51) and (52) we 
see that 



(53) 



vlit,0) = 0, z;i(t,0)= 7172^3=" on(0,r). 



Consequently, ( 40 ) reduces to 



(54) 



f vl=vMt,x)eiO,T)x{0,l), 
v^ = vl,,{t,x)e{0,T)x{0,l), 
v^{t,l)=vl{t,0) = 0,te{0,T), 

wi(i,o) =7172^3°°, ie(o,r), 

v^{0,x) = W°°'\x), 1)2(0, x) = W°°''^{x), x& (0,1), 
[ xi{t) = X2{t) = 0, x^{t) = Zf, t € (0,T). 



+ 00 



We recall that if W^'^{x) = J2 afcCos(i^^^^) and W^'^x) ^ ^ 6fe cos(^^^^^), then 

k=l 

v^ can be expressed , for any (i,x) £ (0,r) x (0, 1) as 



+00 

fc=l 



V {t, x) — y. '^k COS 
fc=l 



{2k+l)TTt 



2bk f (2k + l)TTt\ 1 f(2k + l)nx 

-— ^ sm cos 

(2fc+l)7r V 2 J\ \ 2 



At a; = 0, we have 



OO 

k=l 



ak cos 



/(2fc + l)7rt 
V 2 



2bk /(2k + l)TTt 

■ sm ' 



(2fc + l)7r V 



te (o,r). 



10 



Multiplying this last equality by cos ( - — ^ ) j f^^' ^W fc G N, and integrating on (0, 2) we 
obtain 

Ofc = 0, fc e N, 

and thus, 

(55) W^°°'i=Oon(0,l). 
In particular, 

(56) W^°°'i(0)= 7172^3^=0 
and 

«,i(i,0) = f;6,cos((^^^±^),ie(0,T). 

k=l 

Parseval's identity then implies 

.2 oo 



(57) 


/ \v}(tM'dt > 
Jo 


•^0 fe=l 


Thus, 


from (54|, we finally get 


^oo,i ^ ^oo,2 ^ on (0, 1) 


and 




Z^ = Z^ = Z^ = 0. 



This is a contradiction with ( 36 ) and finishes the proof of Theorem [SJ 

4 Numerical simulations 

In this section, we illustrate the efficiency of the above source estimation algorithm through 
numerical simulations. We consider the system (IT]) with source term q := x — x^, together with 



the estimation algorithm (26l-(27) with initial estimate q = 0. We fix the observation horizon to 



T = 3 and consider 50 iterations of the estimator (|26[)-(27|. The simulations of Figures land 2 
illustrate the performance when we have added 10% white noise on the measurement output and 
where the observer gains 71 and 72 are chosen to be 

71 = 1, 72 = 1/2. 

The numerical simulations have been done through a finite difference method where the time and 
the space are discretized simultaneously. We have chosen a spatial discretization with 20 steps 
(Ax = .05) and a CFL coeflricient of .005 (At = 2.5e - 04). 

5 The A^-dimensional case 

The aim of this section is to show that the Claim [l] can be easily extended to the iV-dimensional 
case. More precisely, let fl he a bounded domain in M", with smooth boundary F. We assume 
that F is divided in two parts F = Fq U Fi, with Fq n Fi =0 and Fq satisfying the geometrical 
optics condition [5]. That is, there exists a positive constant L-pg such that every generalized 
geodesic (the rays of geometrical optics) of length greater than Lpo passes through Fq at a 
nondiffractive point. Finally, let T > 0, w e E and q G H^{n) n Hq{Q). We denote W := {w £ 
H^{Q) s.t. w\ri — 0}. We consider the following system 

Vtt - Av ^ q{x) cos(cjt), {t, x) e (0, T) x 17, 

v{t,x) =0, {t,x) e (0,T) X F, 

(58) <( vlo,x)=vt{0,x)=0,x(zn, 

dv 
y(i,x) = — (i,x), (t,a;)G(0,r)xFo, 
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Iteration number 







Spatial coordinate x 



Figure 1: The source estimator q is traced after each iteration; as it can be seen, the estimator 
after 50 iterations has converged towards q. 
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Figure 2: The first plot illustrates the decrease of the L'^ distance between q and q after each 
iteration and its convergence to zeros. The second plot illustrates the decrease of the Lyapunov 



function defined in ( 46 1 
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where {u, ut) represents the state of the system, y is the output and we set 



"(i,a:) = Y^iy^{^x)j^{t,x), {t,x) e (o,r) X r, 



On 



dxi 



where i'{x) — (i^i(a;), ..., Vr{x)) is the outward unit normal to F at x. 
Then we have the following result 

Claim 9 We can construct efficient observers for which the back and forth algorithm, is conver- 
gent and which allow to reconstruct the unknown source term q{x), using the measurement output 
yit) over the time-interval (0,T), where T > Tabs 'is any time strictly greater than the minimal 
observability time of the wave equation observed at Fp /^/. 

To prove this claim, we proceed as in the proof of Theorem [Tl We begin by introducing an 
equivalent system with w — uj'^v-\-Vtt, so that the unknown parameter is transferred to the initial 
state and we can easily see that the output of this second system, Y, is equivalent to the output 



of Q in the space H^iO,T;L'^{ro)). As before, we extend (w^ , w'^ , z^ , z^) to ((0,+oo) x fl)^ x 
((0,+oo) X Fq)^ and we propose exactly the same observer as in 1-dimensional case. The error 
dynamics are therefore given by 

' Wl = W^, {t, x) G (2fcT, (2fc + 1)T) X n, 
Wi = AW\ {t, x) G (2fcT, (2fc + 1)T) X n, 
(59) \ Wl^~WMt,x)€{{2k + l)T,{2k + 2)T)xn, 

Wi = -AW\ (t, x) G {{2k + 1)T, {2k + 2)T) x Q., 
W^'(t,a;) = 7iJi(i,3;)+7i72^^(i,x), {t,x) G ((2fc + l)r, (2fe + 2)r) x Po, 



and 



(60) 



' Zl{t,x) 


Z?{t,x) 


Zl{t,x) 


ZUt,x) 


ZUt,x) 



^Z\t,x)--i2Z\t,x), {t,x) G (2fcr, (2fc + l)r) xFo, 



— —u) Z {t, X 



dn 



ft, x), (t, x) G {2kT, {2k + 1)T) X Fo, 



-Z^{t,x)-'y2Z\t,x), {t,x) G ((2fe + l)r,(2fc + 2)r) x Fo, 



df: 



-{t, x), {t, x) G {{2k + l)r, {2k + 2)T) X Fo, 



Z'-{t,x), {t,x) G ((2fc + l)r,(2fc + 2)T) xFo, 
Here, the unique solution is defined on the state space 

{WfW^) e C([0, +oo);H^{n) n W) X C([0, +oo); W), 

{ZfZ^,Z^) e ili([0,+(X));L2(Fo)) x L^{[0,+^);L^{To)) x H^[0,+^); L^To)) 
Furthermore, we have the energy estimates 



(61) 

|VVF'(i,.)lL2(n) + ilF'(t,.)li2(n)+7i|^'(t,.)l^(r„)+Wl^'(i,.)li2(r„)+27i72^'l^'li2((o.i)xro) 

|VtFoii2(n) + |W^o 1^2(0) +7i|^olL2(ro) + 7i^^l^olL2(ro) 
and 



(62) -(i W/(t, .)il2(o) + \VW'{t, .)|i2(„) + 7i^1^'(t, ■)nLHr„) 
1 



+ r^' dn 



dWx , ,,2 ,1 2 1 ,^2 1 2 , 2(^2/, M 2 

{'t,x)\L2(ro)+ 2'^n2^ \Z lL2((o,t)xrn) +711^ \Z {t, ■)\l2(t^,) 



< C„(|Wo Ih2(s2) + Wo\H^(n) + l^olL2(ro) + l-^ol )l2{To)^ 
where C^ denotes a positive constant which only depends on uj. 
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Finally, assuming T longer than the observability time of the wave equation, we have the 
following result. 



lim {W'^{2nT, .), W^{2nT, .)) = (0, 0) in W x L^{VL), 
lim (Zi(2nT,.),^^(2nr,.),^3(2nr,.)) = (0,0,0) in ^^(ro) 

n—¥-\-oo 



(63) 
(64) 

In order to prove (|63]) and (64), one may proceed exactly as in the proof of Theorem Is] The 
main difference is in the fact that by the energy inequalities, we have Z^ S L°°(0, +00; Z. (Fq)), 
and thus (up to a subsequence of (0(n))„>o), there exists Z]^ € L^(Fo) such that 



(65) 
when n 



and 



Z^{2(t){n)T, .) -^ Zl^ weakly in L^{Vo) 



\Zl\mv.-^ < hminf |Zi(20(n)r,.)U2(r„ 



We aim to prove that Z^{(t){n)T, .) converges strongly to in i^(Fo) We start with introducing 
the following Lyapunov function 

V{W\W\Z\Z^){t):=\{[ {\VW\t,x)\^ + \W\t,x)\^)dx 

(66) ^ ■''^ r ^ r ^ 

+ wM \Z^{t,x)\^dx+ I \Z'^(t,x)\^dx),te{Q,T). 



Indeed, V{t) > 0, t > and applying (59 1 and (60), one can compute, for any time t G (0,r), 



(67) 



dt 



V{W\W^,Z\Z'^){t) = -Lj'^ / \Z\t,x)\'^dx<0 



The function V is positive, decreasing. Consequently there exists I > such that 



(68) 



lim V{W\W^,Z\Z^)(t) = /. 



From (67) and (68), 



(69) 



+00 



\Z\t,.)\l,^^^)dt = l-V{Q)<+^. 



The idea is to apply the Barbalat's lemma. In this aim, we need to prove that -^\Z^{t, Oliajr ) 
is bounded on R+ (which implies the absolute continuity). This can be done through simple 



computations based on ([59| and (60). Consequently, from Barbalat's lemma, 

(70) lim Z\t,.) =OmL^{ro). 

t—>--\-oo 

In the same way, there exists Z^ G L^(ro) such that 

(71) Z'{2^{n)T,.)-^ZlmL\ro) 

when n — > +00. Once again, we want to prove that Z^ is indeed and furthermore that, the 
convergence is actually strong. We first remark that f{t) := |Z^|^2(p > is absolutely continuous 
on IR-|_. Indeed, for any t E M+, 



(72) 



fit) = -2u' 



Z^Z^ dTn - 2 



xodW^ 



dn 



dTo, 



and therefore, (61) and (62) imply the boundedness of f'(t) on M+. 
Now, one has, at leastTormally, 

f{f z^z^dVo) = f z\-io^z'-'^)dro+ f iz^-j2Z^)z^dro 



IZ'l'dTo 



z'^dro+ f \z'?dv,~^J z' 



Z'^dTo. 
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In other words, 






= jif z'z^dro) + oj^ f l^'l'dro + 72 / z'z^dVo 

^Fq ^Fq -/Fq 



which finaUy gives 



g{t) + h{t) < ^(f Z'Z'dTo) + {u' + 4) f \Z'\^dTo. 



'7iaw^ 



Here, g{t) := — jp Z^^^dFo which tends to when t —> +00 from (62 1 (which impHes the 
boundedness in i'^(ro) of ^^^ ) and (70). Also, h{t) :— ^ Jp \Z^\'^ dFo^ is positive and absolutely 



continuous from ([72|. Integrating the last inequality on M+ and applying (61) (which implies 
the boundedness in L^(ro) of Z^), (701 and (69) we obtain. 



+00 



{g{t) + h{t))dt < 



-^2 r+00 

Z\0, x)Z^{0, x) dVo + {lo^ + ■^) 

Fn ^ Jo JTa 



IZ^fdTodt 



< +CX3. 

Barbalat's lemma implies that 



h{t)^- iZ'^l'^dTo^O ast-^+00. 

In particular, we have 

Z^{2(j){n)T, .) ^ in L^(To). 

Now, from the boundedness of \Z^{2(f){n)T, .)|i2(ro) and \W^{2(j){n)T, .)|L2(ro) (which is a conse- 
quence of the convergence of W^{24i{n)T, .) in W), and the fifth equation in ( [59| , we also have, 
up to a subsequence, that there exists Z^ e iJ^^(ro) such that 



(73) 



hin Z'{2(f>in)T, .) = Z^ in H-\To). 

n— J-+00 



As in the proof of Theorem [sj we now define {v\,v^) G C{[0,T]-W) x C{[0,T];L^{n)) and 
(C^C^C^) eH\0,T;L^iTQ)) X L^{0,T;L^{Ta)) x H^{0,T; H-\To)) as the solution of 

^ v} = v^, (i,x) e (o,r) X n, 
v^ ^Av\ it,x)e (o,T) xn, 

v' - liC + 7172C', {t, x) G (0, T) X To, 
Ct'-C'-72C\(i,a;)€(0,T)xro, 

C? = -^V(i) + |;^, (i,x)e(0,r)xro, 

C3 = C\(i,a;)e(o,T)xro, 
v\o,x) = M^i(a:), «2(o,x) = wlix), xen, 
{ eiO,x) = 0, C'(0,x) = 0, C'(0) = Zi(x), X e To. 

By continuity of flow with respect to the initial state, we have 



(74) 



(75) 

lim {W\2(j){n)T + t, .),W\2(f>{n)T + t, .)) ^ {v\v^) in L°°((0, T); W) x L'^{{0,T); L^{n)), 

n— J- + 00 

(76) 



lim {Z\2^in)T + i, .), Z2(2</)(n)r + i, .)) = (C\ C) in i°°(0, T; L2(ro)), 



n— »+oo 



(77) 



lim Z-^(20(7i)T + t, .) = C-^ iiiL°°{0,T;H-\Ta). 

n— )-+oo 



Following the same strategy as in the proof of Theorem [Sl we have 



C^^o, c'^o, cHt,x)^zlix), 



dn 



on (0, T) X To 
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Consequently, we obtain 



(78) 



vl = w2, {t,x) e (o,r) X n, 

vf = Av\ it,x)e (0,T) xn, 

v^{t,x) = 0, {t,x) e (o,r) X Ti, 

— {t,x)^0,{t,x)e(Q,T)xro, 
v'\t,x)^^i-l2ZUx), {t,x) e (o,r) X To, 

I v\o,x) = wl{x), vHq,x) = wl{x), xen, 

whose solutions are defined on the space C([0,T];W) x C([0,T]; L^(ri)). As we have assumed 
T > Tobsj we can find a sufficiently small t such that T — r is also greater than Tobs- T ^t > Tobs- 
We consider now the function 

(■0^ (i, a;), -0^(i, x)) := (w^ (t + r, x) - w^ (i, x),v'^(t + t,x)- v^ (i, x)) 

defined on the space C([0, T - t]; W) x C([0,T - r]; ^^(f^)). The functions (V-H^^a;), V'^(i, x)) 
satisfy the following equation: 



(79) 



C ^1 =02^ (t,a;) e (0,T-t) X 17, 
i)1 = AV/ , (t, x) e (0, T - r) X f], 
V'^(t,a;) =0, (t,x) G (0,T-r) x T, 

^(i,x) = o, (i,x)e(o,r)xro. 



I V'(0,x) = «1(t) - lyi(x), 02(0, a:) = v^{t) - W^^(x), x e fi, 
But as T — T > Tobs this necessarily means that 

V'^ = V^ = 0. 

As this is through for any r such that r e [0,Tobs — T), we necessarily have that v^ and u^, 



solutions of (781, are constants of time: 

v\t,x) = WUx), v'{t,x) = Wl{x). 



Inserting this inside ( 78 ) , W^ (x) = and W^ (x) must satisfy the following elliptic equation: 



(80) 



AW^ = 0, x G f7, 

W^{x)^0,xeTi, 

W^{x)^ji^2ZUx),x£To, 

-^{x) = o,xero, 



The unique possible solution to this elliptic equation is trivially W^ = and Z^ = 0. This ends 
the proof of the convergence and therefore implies the Claim [9] ■ 



6 Appendix 

This appendix is devoted to the proof of Proposition [7j 

6.1 Existence part 

We first concentrate on the existence of a solution to the forward system. More precisely, our 
aim is to prove the following result 
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Proposition 10 There exists Tq > such that for any 71, 72 > 0, for any {q'^,q^) S H'^{0, 1) H 
H^{0,l)xH^{0,l) and{^°,^^,c()£R^, there exists {v^ , v^) G C([0, Tq]; ^^(0, l))xC([0, Tq]; ff-i(0, 1)) 
and (Ci,6,6) G H^{0,To) x L2(0,ro) x H^{0,Tq) solution to 



«2, (i,x)e(o,ro)x(o,i), 



(81) 



t;,^=«i,, (t,a:)e(0,ro)x(0,l), 

z;i(i,l) = 0, z;i(t,0) = 7i6(i) + 71726 W, t e (0,To), 

vi(0,x) =q°(x), v^{0,x) ^q^{x), xG (0,1), 



(82) 



Lit) = Ut) - 72Ci{t), t ^ {0,To), 
Ut) = -uj^Ciit) + vl{t,0),teiO,To), 

6(t) = a(i),ie(0,ro), 
[ 6(0) = e?, 6(0) = CS, 6(0)^^3"- 



Proof of Proposition [To| Let 71, 72 > 0, (g°,g^) G ii"^(0, 1) n i?o(0,l) x H^{0,1) and 
(^1:^2,^3) e M^ It follows from [3 Theoreme 4.2 page 46] that for any T > 0, for any / e 
H^(0,T), there exists a unique solution (?;^,w^) to 



(83) 



vl^v^,it,x)€iO,T)x{0,l), 
«2=.^;i^, (i,a:)e(0,r)x(0,l), 
i;i(t,l) = 0, «i(t,0) = /(i), te(0,T), 
v^{0,x) ^q"{x), 1)2(0, x) =g^(x), xe (0,1), 



with {v\v^) e C([0,r];L2(o,l)) X C{[0,T];H-\0,1)). For any T > we define J"t on H\0,T) 

by 



J-T : / e i?i(0, T) K^ 716 + 71726, 



6(0 = 6(i)-726(i), te(0,T), 
6(i) = -^'6(0 + «i(^,o),ie(o,r), 

6(i)=6(0,<G(0,T), 

6(0) = e?, 6(0) = 6°, 6(0) = 6°, 



(84) 

where 6 and 6 satisfy 

(85) 



where v^ € C([0,T]; £^(0, 1)) denotes the solution of (|83|. Our aim, here, is to prove that there 
exists some Tq > such that Tto has a fixed point. 
We first state and prove the following Lemma 



Lemma 11 The solution v^ to (83) moreover satisfies the following regularity property 
(86) \vl{t, 0)|i.(o,^) < 2(4r2 + 3)|/|2,,(o J,) + 2(2 + T){\q%^^,,^ + \q'\h^o,i))- 
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Let t £ (0,T). We multiply the second equation in ([83| by w^ and 



Proof of Lemma 

integrate on (0,t) x (0, 1). An integration by parts gives 



k'(^,-)li2(o,i) + l«'(^,-)li2(o,i) 



«2(s,0)f^(s,0)ds + |ql|i2(o,i) + |<?°|i2(o^i) 
< eil«'(^,0)|i.(o^j,) + ^|«i(t,0)|2,(Q^^) 

+ \l Il2(0,1) + I9xIl2(o4), 

where ei denotes a positive constant that will be chosen later. This last inequlity holding for any 
time i € (0, r), we finally obtain 

sup Qv^it, .)|i.(o 1) + \vlit, .)|i.(o,i)) < eik'(i, 0)1^2(0.1.) 
(87) ^ 

+ -\vlit,o)\hio,T) + \<i'\h{o,i) + kSli2(o,i). 
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Now, we multiply the second equation in (83) by {x~l)vl and integrate on (0,T) x (0, 1). Several 
integrations by parts give 

T' 1 T"' T^ 1 

hf f \v^{t,x)\^dxdt~ f \v^{t,0)\^dt+ f f \vl{t,x)\^dxdt) 



Jo 



Jo 



+ {x~l)vl{T,x)vl{T,x)dx- {x-l)vl{0,x)vl{0,x)dx- - \vl{t,0)\^ dt ^ 



Using the boundary conditions in ( 83 1 this leads 

\vl{t,0)\hio,T) < f f {\vHt,x)fdxdt+\vl{t,x)\^)dxdt-\f\l.(^o. 



^0 



T) 



+ \vliT, .)li^(o,i) + l"'(^' ■)\h(OA) + \q°\Ihoa) + k'li^(o,i)- 



Combining now (87) and (|88|) and taking ei := AT, we obtain 



+ (l + r)(k2li^(oa) + k'li^(o,i))- 



We use one more time ( |87[ ) with ei := 4 and we finally get 

l^lit, 0)fmo.T) < 2(4r2 + 3)|/|^,(o,T) + 2(2 + T)(|gO|i.(o,i) + |g^li2(o,i)), 
which ends the proof of Lemma [TT] ■ 

Now, we prove the following proposition ensuring that the map Tt is well-defined. 

Proposition 12 For any T > 0, Tt '■ H^{0,T) — )• H^{0,T) is well-defined. More precisely, for 
any f G H\0,T) there exists a unique solution (^1,6,6) e H^{0,T) x L^{0,T) x H^{0,T) to 



(85) 



Proof of Proposition [l2| Let T > and / G H^{0,T). We define yl e 7W2(M) and B e 



^:= 



-72 1 

-LJ^ 



B := 







where {v ,v ) denotes the solution of (83) associated to /. It is clear that (Ci,'?2,'^3) defined by. 



(89) 



a(i) 

6(i) 



-,<>! 



($)+l ^^'~'^^B{s) ds and 6(t) :- ^3° + ^ ^(s) ds 



is the unique solution of (l85|. Let us prove that (^1,6,6) e H^{0,T) x L'^{0,T) x H^{0,T). 



Simple computations lead to 



l(ei,6)li: 



(90) 



(0,T) 



*^' ^° Mi^(o,T) + l/ e(*-^)^i3(.)dsU.(o,T) 



< le I to 







< CTe^^(|e?P + |C."P) + C^ te^*^*k(.,( 



l^dsl dt 



where we recall that C = C(a;,7i,72) denotes various positive constants which may vary from 
line to line and only depend on w, 71 and 72. Applying Lemma ITTl (90) becomes 

|(^i,6)li.(o,T) < CTe^''m\' + \e2\') + CT'e^^\vlit,0)\l.^„^^^ 
(91) < CTe^^(|^?p + |^2°n 

+ 2Ce^^T^ ((4r2 + 3)|/p^,(„^j,) + (2 + T)i\q%,^,,^ + \q'\h(^o,i)) 

From the last inequality, (^1,^2) G i^(0,T). Consequently, from the first line in (85), 
(92) ^ieH\0,T). 



The fact that ^3 e H^{0,T) finally comes from (89) and (92). This ends the proof of Proposition 

nn ■ 
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Now, we prove the following result, implying the existence part of the Proposition [7] 

Proposition 13 There exists sufficiently small Tq only depending on uj, 71 and 72, such that 
Ttq admits a fixed point in H^{0, Tq). 



Proof of Proposition 13, For any i? > we call Br := {/ G iJ^(0,T) s. t. |/|ffi(o,r) < R}- 

In order to prove Proposition |13[ we are going to prove the two following points, for sufficiently 

small To, 

i) there exists i?o > such that J^To{Bro) C -Bflo, 

ii) J-jb is contractant. 

Point i): 

Let i?o > be such that 

(93) Ro > \e,\' + IC2P + legP + \<}%Ho,i) + \l%i(o,iy 



and let f e Br^. From (85), M and (91 1, 



I^3|h1(o,T) ^ I^3|l2(o,T) + I^iIl2(0,T) 

< T\C!\'+T'\Ci\h(o,T) + \^i\h(o.T) 

< T\e,\' + C{l+T^){Te^^i\e^\' + \e2\') 

+ 2Ce^^T2 ((4r2 + 3)|/|?,,(o,^) + (2 + T){\q'^\l.^Q^,^ + \q'\h^„^,^)) 



(94) 



By straightforward computations, there exists sufficiently small Tq, only depending on (w, 71, 72), 
such that 



I716 +7l 7261^1(0, To) < -Ro, 



which implies the Point i). 
Point ii): 



Let /i, /2 e Br„. We denote by {v^,v^) and (Ci,6,6) (respectively by {ri^,r]'^) and ((1,(2,(3)) 
the solutions to (83) and (85) associated with /i (respectively with /2). Then {v^,v'^) := {v^ — 
7f,v^-ri^) and ((1,6,6) := ((1,(2,(3)- ((1,(2,(3) satisfy 



(95) 



and 



~2 _ ~1 



^i^(t,a;)G(0,T)x(0,l), 
vl,,{t,x)e{0,T)x{0,l), 



vHf 1) = 0, z)i(t,0) = /i(t) - /2(t), i e (o,r), 
v^{Q,x) = 0, w2(o,x) = 0, a;e (0,1), 



(96) 



(1 (i) -(2(0- 72(1 (i),te(0,T), 
(2(i) = -c^2^;(t) + «i(t,o),ie(o,r), 
(3(i)=(i(i), <e(0,T), 
I (1 (0)^0,(2(0) -0,(3(0)^0. 



Adapting Lemma 11 (91) and (|94[), one easily obtains 



|7i(i+7i72(3ll,i(o,T) < Ce^^T2(4T2 + 3)|/i-/2|?,,(o^^), 

where C is a positive constant only depending on (01,71,72). Therefore for sufficiently small Tq 

only depending on (01,71,72), Ce^^°TQ{4TQ + 3) < 1 and the Point ii is proved. 

From Point i) and Point ii), J^Tq admits a fixed point in _ff^(0,ro): this ends the proof of 



Proposition 13 We emphasize the fact that the time Tq defined in this proposition, does not 



depend on the initial conditions (Q°,g^) € ^^(0, 1) n H^{0, 1) x H^{0, 1) and ((?,(^,(§) e 
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In the same way, one obtains the following result (existence of a solution to the system 
backward 

Proposition 14 There exists Ti > such that for any 71, 72 > 0, for any {q'^,q^) € H^{0, 1) H 
H^{0,l)xH^{0,l) and{^1,^^,^f) eR^, there exists {v^ , v^) £ C([0, Ti]; ^^(o, l))xC([0, Ti]; iJ-i(0, 1)) 
and (6,?2,6) e H^iO,Ti) x L^{0,Ti) x H^{0,Ti) solution to 



(97) 



vj ^-v^, it,x)e (0,Ti) X (0,1), 

vf^-vl,,it,x)e{0,T,)x{0,l), 

v\t,l) = 0, v\t,0) = 716 (i) +71726 W, i e (0,Ti), 

t)i(0,a;) = (70(a;), ^2(0,0;) = q\x), x& (0,1), 



(98) 



CiW = -6(0-72ei(0,te(0,Ti), 

e2(<)=c^'ei(i)-^'^(i,0), iG(0,Ti), 
6(t)=6(i),ie(0,Ti), 

I ei(o) = e?,e2(o) = e2°,6(o) = e3°- 



We are now going to prove that times Tq and Ti defined in Proposition [TO] and [14] can actually 
be chosen arbitrarily. 



6.2 Global well-posedness 

In the context of an established local well-posedness result, it suffices to prove the following global 



a priori estimate for smooth solutions to system (81 1-([82| (the fact that in this case it will also 



hold for the solution of system (97)- (98 1 is obvious) 



1^1(1^0,. )li.(o,i) + l^'(T^o,.)li.(o,i)+7i|6(ro)P+7i^'iei(ro)P + 27i72u;2|Cili.(o,To) = 



and 



lilvtin, .)li^(o,i) + l^'(^o, .)li=(o,i) + 7i'^'iei(To)n + Ij^ivHtoM' 

To 



(100) 



where Ci^ 



1 2 

+ 7T7i72W 



Ms)\'ds + jiL,'MTo)\'< 



Cw,7l,72(k lff2(o_i) + |g |//1(04) + I6I +1^21 ): 

denotes a positive constant which only depends on w, 71 and 72. One can note 



that these inequalities also ensure more regularity for solutions to (81 1- (82 1 



Step 1: Proof of (99 1 



Let To be defined as in the Proposition [lO] and let (w^w^Cl:6,6) € C([0, Tq]; ^^(o, 1)) x 
C([0, To];H-'^{0, 1)) xi/i(0, To)xL^{0, To)xH^{0, Tq) denote a solution to (|8lJ. Wc first multiply 
the second equation in (81 1 by w^ and integrate on (0, 1). A straightforward integration by parts 
gives 



1 



dt{ \v'\' + \vin+vl{t,Oy{t,0)^0. 



This equality implies, from the boundary conditions in (81) 
»i 



1 



dt{ / y\' + K\') + (6 + ^'a)(7iei + 7i72?3) = 0, 



hence 



1 



dt{ (k'P + |t'^nd:E + 7i|6l'+7iw'l6n+7i72^'iei|' = 
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and (991 follows. (We first obtain the previous equality in the case where (w^,w^) is assumed to 
be smooth enough and we then use a density argument to obtain it in the general case). 



Step 2: Proof of (100) 



Now we derive the second equation in (81 1 with respect to the time variable and we integrate 
on (0, 1). (In the same manner, we first assume that (w^, w^) is smooth enough and we then use 
a density argument to obtain it in the general case.) It follows, 



(101) 



^d,J^ i\vt\' + \vi\') + vi{t,0)vtit,0)^0. 



Using the boundary conditions in (81 1 we compute 



vUt,0)vHt,0) = vlit,0h,^2 



= 7i«^(i,0)(-^2^i+«i(t,0)) 



We integrate ( 101 1 on (0, To) and use the initial conditions in (81 ) and the last identity to obtain 



2(l9LlL(o,i) + k^li^(o,i) +7i|'Z°(0)n. 



Using one more time the boundary conditions in (81 1, several integrations by parts give 



-71C.2 f \^{t)vl{t,Q)dt - 7^ / '' i^{t)vl{t,Q)dt^-i^uj^[(,^{t)vl{tMl° 



(103) 



Now, 



7^ / Um2{t)+uJ^Ut))dt'li^^MtX{tMl 
Jo 



To 



ll^ 



{\^i{Tor-\c'in-ii^' a(i)6(t)dt 



To 



+ii^^[iim2{t)]l"-ii^'Mt)vi{t,o)]o"- 



-71^2 f ° ^■^(^t)Ut) dt 
Jo 



(104) 



-716^2 / (^^(^) _ ^^^^(^) _^ 72ei(i))6(0 dt 

T 

-^{MTo)\'-\e2\')+lil2co'J^"Mt)\'dt 

-iiiW I \i{t)i2{t) dt. 

Jo 



Combining now (102|-(104), we obtain 

2(K^(ro,.)li^(o,i) + k'(7^o,.)li^(o,i)+7ik^(To,0)p + 7io.^lCi(To)P+7ic^'l^2°n 

+ 7172^' I ° 16(01' rfi + 7i^'l6(ro)p-7i72w'Ci(ro)6(ro) + 7iw'e?92(o) = 

(105) -^0 1 

2(kLli^(o,i) + k^li^(o,i)+7i|'Z°(0)|'+7i^'lC?l' + 7i^'l6(To)n 

+ 7i7l^' / \i{t)Ut)dt + uj\'i\'~lil2U^^eiil + li^^UTo)vl{nM- 
Jo 
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Using the following estimates 



7i72^'a(To)6(To) < ^^^(lei(To)P 



7i-^C?9°(0) < ^(|C?I 



92(0)n, 



2 2 

7i72^ 



eiW6Wrfi<^^klli2(0,To) 



71 72^ 



1^:'^ 



2Il2(0,To)' 



7l72C^'?i'^2" < 



-mi' + iein 



^iUj^^i{To)vl{To,0) < 2j,uj^\^i{To)\' 



71 



^^(^0,0)1^ 



together with (99) and (105), we finally get the existence of a positive constant C^j^-yj^^j which 
only depends on cj, 71 and 72 (we underline the fact that 6*1^,71,72 does not depend on the time 
To) such that (105) becomes 



2(k'(To,.)li^(o,i) + l^'(7^o,.)li=(o,i) +7ic^'ICi(To)n + ^k(ro,o)p 



+ 







\ut)rdt+jiuj'\uTor< 

CLj,7i,72(k l//2(o,i) + k lffi(0,l) + l^il +IC2I ) 



This ends the proof of the energy estimates. An immediate consequence of these energy estimates 
is the uniqueness of the solution (v\w^6,6,C3) e C{[0,To]; L^{0,1)) x C([0,To]; if-i(0, 1)) x 
H^{0,To) X L2(o,To) X fl^2(-o^j|j) to ([|l])-([^. Another consequence is that this solution is 
actually more regular, more precisely, 

(i;\t-') € C([o,ro];ij2(o,i)ni/,i(o,i)) X c([o,ro];^.'(o,i)). 



Finally, the fact that the right members of both estimates (99) and ( |100 ) do not depend on 
the time Tg ensures that the solution to (81 ) does not blow up at time t = Tq and can also be 



extended at the right of this time. In other words, the solution to (81 ) exists, for any time T > 



In the same way, we obtain the same result for the system backward (system (97)-(98)). 
Existence and uniqueness of a solution (w\ w^,CiiC2,C3) € L^^((0, +00); _ff^(0, 1) n H^{0, 1)) x 
^l?c((0>+oo);ff.^(0,l)) X Hi„^(0,+oo)) X Ll^{0,+^)) X Hl^{0,+cx^)) to system ^-^ im- 
mediately follows. 



Let us now prove that energy inequalities are global in time, i.e., let us now prove (32) and (33) 



6.3 Proof of energy inequalities 



Step 1: proof of (32) 



Let t G (0, +00). We assume that there exists fco > such that 2koT <t< (2fco + 1)T. (The 



proof of (32 ) in the case were (2fco + ^)T < t < (2kQ + 2)T readily follows and we omit it.) One 
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can compute easily, using (30) and (31), 



fcn-l 



f.(2fc+l)T fl 



/ / v^v^dxdt = Y.^ / ^lydxdt^ / viydx 

Jo Jo .,_n J^kT Jo J(2k+l)T Jo 



k=0 



+ 



•.dt) 



vl^v'^ dx dt 



l2koT Jo 
fcn-l ,.(2fc+l)T ^(2fc+2)T ^1 

w^j,tij^ dxdt + / wiJ;2,Wt^ dx di 

l2kT J{2k+1)T Jo 



/c=0 



+ / / vlxvldxdt 

2koT Jo 
t rl 



fco-1 ^(2fc+l)T 



-1 pt pL r^u-i p(iK+L)l 

0/ dti\vl\^)dxdt-J2 (6+^.2^1)716^* 

^ Jo Jo j.^Q J2fcT 

fco-1 ^(2fe+2)T <.t 

V / (-6 + c.'a)(-7i6) d<) - / (6 + c^'6)7i6 dt 

k=0 "'(2fc+l)T J2feoT 

t fl 

"'0 

fco-l ^(2fe+l)T 

+ 12 I (-7166 -7iw'Ci(ei+ 726)) rfi 



fc=0 



2kT 



fco-1 i-{2k+2)T 

y] / (-7166 + 7i'^^Ci(-6 - 726)) dt 

fe=0 •''(2fc-+l)T 

+ / (-7166 - 7i^^Ci(Ci + 726)) dt. 

J2koT 



The end of the proof of ( 32 ) follows easily. We now focus on ( 33 ) 
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Step 2: proof of (33) 



We still assume that we are in the case where 2koT < t < (2fco + 1)T- Then, 

t /.I fco-J- /.(2fe+l)T /.I /.(2fc+2)T fl 



"'0 



v?tV? dx ds 



fco-l f(2k+l)T i-l 



fc=0 



2kT 



v^^iVi dx ds 



(2fc+l)T Jo 



Vxxt'Vt dx ds) 



vl.^f.Vf dx ds 



2feoT ^0 
t rl 



ko-l i-{2k + l)T 



2feoT 



-J dtilvl^ndxds-Y, vUs,OK{s,0)ds 

^ Jo Jo ^^Q J2kT 

ko-l i-{2k+2)T j-t 

Y, / vUs,0)vUs,0)ds- I vUs,0)vUs,0)ds 

l^^Q J(2k+1)T 

f dt{\vl\^)dxds 
^ Jo Jo 

fco-1 „(2fe+l)T 

- Y. I vl,{s,0)j,{^uj'^,+vl{s,0))ds 



1 



fc=0 



2feT 



fco-1 <.(2fe+2)T 

^ / z;i,(s,0)7i(-c^'6+«^(s,0))ds 

fc=0 •^(2fc+i)T 
ft 

vUs,Ohi{-w^^,+vl{s,0))ds. 

2koT 



This rewrites 



2 Jo Jo 



71 



2 7o 



id,{\v^\^)+dt{\vin)dxds+^ dti\vl{sM')ds = 7iwM ^;i,(s, 0)6(5) ds 



Now, we still proceed as in the proof of ( 100 ) 



(106) 



710.2 I yl^^s,0)Us)ds - -71^2 y Us)vi{s,0)ds + j,u^[U-sH{s,0)]l 

-71^' E / a(5)(6(s) + ^'ei(s))ds 

fc=0 ^2fcT 
fco-l A2k+2.)T 

-7iw'E / 6(5)(-6(s)+c.2a(s))ds 

fc=0 •''(2'c + l)T 

-7iu;2 /• (,[s){Us)+u:^Us))ds + ^,uj^[Us)vl{sMo 

J2koT 
4 fcp - 1 /.{2k+l)T 

Us)Us)ds 



71 w 



■(i6(t)p-ie?n+7i-^E / 

fe=0 "'2feT 



fco-1 p{2k+2)T pt 

-71^' E / li (s)6 (s) rfs + 71 w' / a (s)6 (s) ds 

fc=0 J(2fc + 1)T J2fcoT 



fc-o-1 



fc-o-1 



-71-^ E Ki(^)a(^)]^f/'^" + 71-^ E Ki(-^)a(^)]S:i!? 



J2feT 
fc=0 fc=0 

.2rt f „\c I „\}t I „, , ,2r 

fco5 



-7iw'Ki(s)6(s)]kT + 7i^'[a(s)«^(s,0)]*. 



fco-l 



fco-1 



LetA:=7ic.2(-^[^i(s)C2(s)](2,'r 



(2fc+l)T 



E Ki(^)a(^)]( w - Ki(5)6(s)]koT)- 



fc=0 



k=0 
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One can compute 



feo-l 



71 "^ 



^ fc=0 



fc=0 



i2kT 



(2fc+2)T 
(2fc+l)T 



Ki(5)-72ei(s)6(5)]LoT 



= 7i^'(-el(i) + ei(o) + 72 E (a((2fc + i)r)6((2fc + i)T)-a(2fcr)C2(2fcr)) 

^ fe=0 

feo-l 

-72 J2 (^i((2fc + 2)T)6((2fc + 2)r) - ei((2fc + l)T)e2((2fc + l)r)) 

fe=0 

+72(6 We2W - 6(2fcor)6(2fcor))) 
= 7i^'(-elW + el(o) 

+72 'E'(6((2fc + l)T)(ei((2fc + 1)T) + 72ei((2fc + 1)T)) - ^,{2kT){U2kT) + i^iii^kT))) 

fc=0 
feo-l 

-72 E (^i((2fc + 2)T)(-ei((2fc + 2)T) ~ 726 ((2A; + 2)T)) 

fc=0 

- ei((2fc + l)T)(-^i((2fc + 1)T) - 72ei((2fc + 1)T))) 
+72(eiW(eiW+72eiW)-ei(2fcor)(ei(2fcoT)+72ei(2fcor)))^ 

= 7i^'(-e|(i) + e2 + 72^1(^)6(0 - 72???2°)- 



Thus, (106) becomes 



71^2 / vl,{s,Q)ii{s)ds 
Jo 



71^ 



fco-l p{2k+l)T 

2 (i6(t)p-ie?n+7ic^^E 

fe=0 

fco-l „(2fe+2)T pt 

• - ' •• 2 



6(s)6(s)ds 



2fcT 



-11^^ y, ei(s)6(s)ds + 7iw' / 6(5)6(5) ds 

j,^g J(2fe + 1)T J2fe-oT 

+7i^'(-el(0 + 6° + 726(t)6(t) - 72C?6°) +71^' [6 (5)^^(5,0)]*. 



Now, for any < k < k 







thus 



7iW 



7ia;- 



-72^ 



^(2fc+l)T ^(2fe+l)T 

7ic^' / 6(5)6(5) ds = 71^' / (6(5) - 726(5) + 726(5))6(5) ds, 

J2kT J2kT 

j.(2k+2)T <.(2fe+l)T 

-7iw' / 6(5)6(5)^5^-71^2/ (-6(5)+ 726(5) +7l6(5))6(s)d5, 

J(2fc+1)T J2fcT 

71^'/ 6(5)6(5) Cis = 71^' / (6(5)- 726(5) +726(5))6(5)d5, 

J2kaT J2knT 



feo-l ^(2fc+l)T fco-1 i-{2k+2)T j-t 

'(E / 'iMUs)ds-Y, 6(5)6(5)^5+/ 6(5)6(5) ds) 

,._n "'2feT .._n J(2k+1)T 



fe=0 



fe=0 



2fcoT 



1 (.t (.t feo-l ^(2fe+l)T 

'(2 i a*(|6(5)Pd.-72 [ |6(5)Pd5 + 7lE 







fe=0 



6(5)(6(5) + 726(5)) rfs 



2feT 



co-l ^(2fc+2)T i-t 

E/ 6(5)(-6(5) - 726(5)) rfs + 72 / 6(5)(6(5) + 726(5)) (i5. 

.._n J(2fc+1)T J2fcoT ^ 
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So, finally. 

, fcp-l j.{2k+l)T fco- 1 i-(2k+2)T /-t 

^ fc=0 "'2fcT j,^(, J(2fe+1)T J2feoT 



k=0 "'2'cT j,^(, J{2k+1)T J2kaT 



Following now the proof of (100), one easily obtains (32). The proof of Proposition [t] is now 
complete. ■ 
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